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Abstract 



o ■ 
O ; 

^ ; Arguments are given favoring possible connection of GRBs with core collapse 
Q ■ of massive Wolf-Rayet stars. We analyze the observed properties of cosmic 
(D ! gamma-ray bursts, Wolf-Rayet fWR) stars and their CO-cores in the end of 
Q ; evolution. WR stars are deprivea of their extended hydrogen envelopes, which 
oo ■ makes it easier for the collapse energy to transform into observed gamma-ray 
^ ! emission. Presently, of ~ 90 well-localized gamma-ray bursts, 21 ones are opti- 
^ ; cally identified and for 16 of them redshifts are measured [z = 0.4-^4.5). The 
> ; observed energy of gamma-ray bursts spans over a wide range from 3 x 10^^ 
^ ; to 2 X 10^^ ergs. There is some evidence that this distribution N{AE) is 

■ bimodal if take into account GRB980425 associated with a peculiar type Ic 
^ : supernova SN1998bw in a nearby galaxy ESO 184-G82 for which AE^ « 10^^ 
O '. ergs. These characteristics of gamma-ray bursts are similar to the distribution 
^ ! of the final masses of CO-cores of WR stars which is also wide and homo- 

geneous: Mco = (1 - 2)M0 ^ (20 - 44)M0. A possible bimodality of the 
^ ; gamma-ray burst energy distribution (£"1 = 10^^ erg; AE2 = 3 x 10^^-!-2 x 10^"^ 

■ erg) is in accord with the bimodal mass distribution of relativistic objects 
^1 [Mns = (1.35 ± O.15)M0; Mbh = {4: ^ 15)Mq). That the supernova 
l; ■ SN1998bw is of the " peculiar Ic" type, atypical for WR collapses (type Ib/c), 

•'-j '. can be related to the rotation of the collapsing CO-core which can make the 
rN ; collapse longer and lead to the formation of a neutron star, the decrease of 

■ the gamma-ray burst energy, and the increase of the fraction of kinetic energy 
transported to the envelope. The expected collapse rate of CO-cores of most 
compact WR stars of type WO in the Galaxy is ~ 10^ per year, which is 
only by one and a half order of magnitude higher than the average gamma-ray 
burst rate in one galaxy (~ 10^ — 10^ per year). Two particular models 
that use WR stars as gamma-ray burst progenitors are considered: the hyper- 
nova model by Paczyhski (1998) and the model of unstable CO-core collapse 
suggested by Gershtein (2000). In both models the allowance of a gamma-ray 
beaming or random outcome of the CO-core collapse due to some instabilities 
permits one to bring the rate of CO-core collapses in accordance with that 
of gamma-ray bursts. We argue that WR stars (most probably, of type WO) 
can be considered as progenitors of cosmic gamma-ray bursts. Two types of 
gamma-ray bursts are predicted in correspondence with the bimodal mass dis- 
tribution of the relativistic objects. Three types of optical afterglows should 
appear depending on which CO-core is collapsing: of a single WR star, of a 
WR star in a WR+0 or a hypothetic WR+(A-M) binary system. In addition, 
we briefly consider a model of gamma-ray bursts as a transient phenomenon 
occurring at early stages of galactic evolution {z > 1), when very massive 
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(M > IOOM0) low-metallicity stars could form. Such massive stars should 
also loose their hydrogen envelopes and become massive WR stars whose col- 
lapses could be accompanied by gamma-ray bursts. WR-galaxies can be most 
probable candidates for gamma-ray burst host galaxies. 

Introduction 

The principal requirements to the central engine of GRBs include (1) The ability to release the 
electromagnetic energy ~ 10^^ ergs during 10-100 s (the typical duration of "long" GRBs, only for 
which X-ray and optical observations are possible; a separate group of short single-peak GRBs is 
much less studied, apart from the fact that they are isotropically distributed over the sky, and well 
may be another phenomenon) and (2) The event rate is on average about one burst per typical galaxy 
(assuming isotropy of the emission and homogeneity of galaxies) per ~ 10'' years. The beaming of 
gamma-ray emission decreases the energy emitted and increases by the same amount the event rate. 

These requirements are met (with different degree of accuracy) by two broad classes of astro- 
physical sources. The first class includes coalescencs of binary NS and/or BH [||]). The fireball 
is generated by neutrino-antineutrino annihilation copiously produced during the coalescence. The 
second class comprises models related to final stages of evolution of massive stars, among which: 

• Collapses of massive stars @, |], 0. An accretion disk is formed around a massive rotating 
black hole during late stages of the core collapse of a massive star; in this model, a narrow 
jet is produced inside the star and punches through the stellar envelope reaching very high 
Lorentz-factors. 

• Electromagnetic model by V.Usov [^], in which the energy comes from the rapid rotation of 
a young neutron star (millisecond pulsar) with a very strong magnetic field. Other models in 
which magnetic field is crucial see in [^. 

• GRB during core collapse of a non- rotating Wolf-Rayet star by Gershtein [Q, in which the 
internal shocks are created due to the collapse non-stationarity and energy is brought away by 
electron-positron plasma. 

At the late stages of evolution, very massive stars lose their hydrogen-rich envelopes (via stellar 
wind or in a binary system) and are observed as helium-rich Wolf-Rayet (WR) stars. They are 
considered as progenitors of lb type supernovae (if strong helium lines are present in the maximum 
brightness) or Ic supernova (if helium lines are weak or absent). If GRBs are directly related to 
collapses to massive stars, there should be some links between properties of WR stars before collapse 
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and GRBs. We suggest that observed broad (and, possibly, bimodal) distribution of GRB energetics 
is related (1) to the broad distribution of final masses of CO cores of observed WR stars before 
collapse and (2) to the observed bimodal distribution of masses of relativistic remnants of stellar 
evolution (neutron stars and black holes). 

GRB association with star formation 

There is a growing evidence that cosmic GRBs are associated with star-forming regions: 

• Optical observations of the identified GRB host galaxies evidence for an enhanced star forma- 
tion rate inside the GRB sites, sometimes by an order of magnitude higher than in our Galaxy 

i- 

• Large column densities 10^^ — 10^^ cm~^ toward GRBs derived from known X-ray and optical 
afterglows of GRBs |jlO|, typical for giant molecular clouds. 



Observed offset distribution of GRBs from their host galaxies statistically suggests association 



with the location of massive stars [|Tl|] 



• Detections of line features in several GRB X-ray afterglows ([|T2| and references therein) imply 
a dense surrounding. 

This favors GRBs association with evolution of massive stars and makes other GRB progenitors 
(such as double neutron star coalescences) less likely. 

Energetics of GRBs with known redshifts 

Table ^ lists GRBs with known energetics. Photometric distances are calculated using a flat Universe 
with = 0.3, = 0.7, Ho = 60 km/s/Mpc. The Table also contains BATSE fluences (50-300 keV) 
and peak luminosities Lp phot/s (from |l20[| ). 

As seen from Table |l| and Fig. |l], ^, the observed GRB energetics spans from ^ 7 x 10^^ ergs 
to ~ 2 X 10^^ ergs. This observational fact is usually explained by a broad luminosity function of 
GRBs (see p3|), although one can construct a self-consistent model with a universal energy 



release of -Eq ~ 5 x 10^^ ergs and a complex beam shape [24]. Adding GRB980425, which is possibly 



associated with pecuhar type Ic supernova SN 1998bw in a nearby ESO 184-G82 {z = 0.0085) |jT9[, 
either evidences for a bimodality of GRB energy distribution or for an extremely broad luminosity 
function (more than 5 orders of magnitude!). 
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Table 1: GRBs with known energetics 



GRB z d], 1028 F^, 10-5 Ref AF, 10^^ Lp, 10" 

cm erg/cm^ erg ph/s/cm^ ph/s 









(10-2000 keV) 






(50-300 keV) 




000926 


2.066: 


5.81 


2.2 f 


13 


3.04 






000418") 


1.118 


2.73 


1.3 


14 


~ 0.6 






0003010^) 


2.03 


5.69 


> 0.05 


15 


~ 0.07 


~ 5 


6.7 


991208^) 


0.706 


1.55 


10 


16 


~ 1.8 






990712 


0.430 


0.85 












990510 


1.619 


4.30 


2.26 [Q 2.0 


8.16 


7.3 


990123 


1.6 


4.25 


26.8 H 23 


16.4 


14 


980703 


0.967 


2.28 


2.26 n 0.75 


2.6 


0.86 


980613^) 


1.096 


2.66 


0.17 [fS 


0.072 


0.63 


0.27 


971214 


3.412 


10.6 


0.944 [Q 3.0 


2.3 


7.4 


970828 


0.958 


2.25 


9.6 [0 3.1 






970508 


0.835 


1.99 


0.317 0.08 


1.2 


0.29 


970228^) 


0.695 


1.52 


~ 0.2 [|T8| 0.034 


3.5 


0.60 


980425 


0.0085 


0.013 


0.32 


19 


7 • 10-6 


0.96 


2.1 • 10-5 



Notes: 

t Flat Universe, Qm = 0.3, Qa = 0.7, Hq = 60 km/s/Mpc 
t Peak fluxes from |2D1 

a) Energy fluence in 25-100 keV; 

b) Peak flux Fp = 3.7 Crab, no fluence published, single-peak proflle, duration 10 s; 
other indirect estimations of total energy in gamma-rays see in [j^T] 

c) Energy fluence for E > 25 keV 

d) Energy fluence for £" > 20 keV 
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Figure 1: Total GRB energy release AE (in 10^^ erg/s) of GRBs with known redshifts from Table ^ 
Figure 2: Peak luminosities Lp (in 10^^ phot/s) of GRBs from Table 

GRB and WR stars 



The arguments favoring the link with WR stars are |25]: 



(1) Energetics The observed GRB energy release AE ~ 10^^ - 10^^ erg is roughly comparable 
with the wide range of CO cores of WR stars before collapse, M[.o ~ 2 - 40Mq, Fig |. The 
total energy released in collapse is Eg ~ GM^/Rc, where Mc in the compact remnant mass, Rc its 
radius. During black hole formation in such collapses without mass ejection the available energy 
can reach 10^^ — 10^^ ergs for the observed mass range of compact remnants, i.e. conversion of 1% 
of the available energy into kinetic energy of shocks with subsequent radiation would be sufficient 
to explain the broad luminosity function. Note that during collapse into black hole without mass 
ejection Rc oc Mc oc Mco and available energy range is proportional to the CO core mass range, 
which is likely to be the case. 



(2) Bimodality of GRB energy distribution and stellar remnant mass distribution. The 

known masses of NS and BH in binaries are shown in Fig. ^. NS masses are groupped ina anarrow 
interval is M^s = (1-35 ± O.15)M0, while BH masses fall in a broader range Mbh = (5 ^ 15)Mq. A 
real gap between NS and BH masses is observed. 

GRBs with low energetics associated with peculiar supernovae type Ic (such as GRB980425) can 
be explained by collapses of bare CO cores of massive stars with significant rotation which causes 
most envelope to be ejected and neutron star to be formed, while collapses of slower rotating cores 
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Figure 3: Masses of final CO cores Mco of WR stars as calclated using observed stellar wind mass 
loss M (X M". Masses of known NS and BH candidates in binary systems are also shown. 

do not accompanied by a significant envelope ejection and lead to black hole formation. In the latter 
case an energetic GRB can be generated with energy proportional to the pre-coUapse core mass. 



(3) Association of GRBs with star-forming regions 



(4) A diversity of the observed afterglows A GRB in a binary system can induce different 
optical phenomena due to illumination of the companion's atmosphere by hard X-ray and gamma- 
radiation [^. This should add some light to the "pure" power-law afterglow from relativistic blast 
wave thus producing a great variety of the observed light curves. These effects can occur in a time 
interval Atopt > D/c after the burst {D is the distance to the optical star from GRB, c is the speed 
of light). Deviations are indeed observed in some bursts (for example, in GRB980326 afterglow three 
weeks after the burst p7|). 

Astronomical observations indicate |23, E3, pOi that about 50% of all WR stars in our Galaxy can 
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Figure 4: Dependence of the masses of NS (circles) and BH (triangles and rectangles) on the 
masses of their companion stars my in close binary systems. Filled circles correspond to radio pulsars, 
filled triangles to BH in X-ray novae. 

be in binaries with 0-star or A-M-star. For example, for WR+0 system V444 Cyg with an orbital 
period of P = 4'^.2 we have D ^ AORQ and the time delay Atopt ~ 100 s, and for parameters of 
WR+0 binary system CV Ser Atopt ~ 300 s. An extremely bright optical emission {V ~ 9™) was 
observed in the famous burst GRB990123 only 50 s after the burst beginning |RT . 



Another example is a peculiar shape of achromatic optical afterglow light curve observed in 
GRBOOOSOlc |]32| , The observed several peaks separated by 2-3 days days can be a manifestation 
of an orbital period in the underlied binary system, for example, through the binary-period shaped 
mass loss before collapse. An alternative explanation by a microlensing event seems less probable. 



Orbital periods of order of several days perfectly fit the observed period range l'^.6-=-2900'^ in WR-l-0 
binary systems (see Table 2 in [pS] ). 
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These arguments favor the GRB-WR stars association, but there is a general requirement which 
should be met by all viable GRB models. The point is that GRB phenomenon should be an extremely 
rare astronomical event. 



Event rate problem 

The GRB event rate per unit comoving volume using BATSE data with Ftr = 0.1 ph/cm^ is 

ncRB ~ lO-^GRB/yr/Mpc^' , 

i.e about 10^^ per year in the average galaxy with a mass of lO^^M©. This is by several orders 
of magnitude lower than the total rate of core collapses associated with SN II and Ibc {TZsNibc ~ 
3 X 10~^/yr/Mpc^, [0]). This discrepancy is usually eliminated by introducing a beaming of gamma- 
ray emission (e.g. [0). It is not excluded that not each SN Ibc is associated with GRB for internal 
reasons. 

The mean formation rate of all types of WR stars in the Galaxy is 

NwR At0 1 _i 
= "^^IV^ ~ 1000^^ • 

i.e. by a factor of 1000 exceeds that of GRBs. The most compact WR stars, so called WO stars, 

are much less frequent (3 out of total 200 are known in the Galaxy), so their formation rate is only 

by one order of magnitude higher than that of GRB. This issue can be solved either by postulating 

generically thin jets or, admitting quasi-spherically symmetric emission, by assuming the existence 

of some "hidden" collapse parameters (rotation, magnetic field, etc.), which was suggested by 

from an independent analysis of black hole formation in binaries. 

In the hypernova scenario by Paczyhski the rarity of GRB phenomenon is also explained by 

requiring an extremely high magnetic field during core collapse of a rotating massive star into a 10 

Mq black hole. 

In contrast, in the model of coalescing neutron star/black hole binaries (which is currently less 
favored by association of all observed GRB hosts with strong star forming regions, see above) the 
event rates varies from ~ 10"^ to ~ 10"*^ per year depending on the binary evolution parameters 



37[] , which is marginally consistent with the observed GRB rate and the event rate problem is not 



very strong. 

GRBs as a transient galactic phenomenon 

There is another possibility to explain the observed association of cosmic GRBs with star-forming 
regions at high redshifts and their extreme rarity. GRBs may represent a transient galactic phe- 
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nomenon occurring at the early stages of galactic evolution, like quasars and AGNs. It is established 
now ||38| that at high redshifts 2; ~ 1 — 2 a violent epoch of star formation in young galaxies occurred. 
It is also known that a lot of cold matter were bound in giant proto-galactic clouds at redshifts 2 > 2, 
which are observed as " Lyman-alpha forest" of absorption lines in quasar spectra. The formation of 
very massive stars 100-500 Mq which final collapse into massive black holes took place at that epoch. 
Such massive star can not form from matter enriched with metals because of pulsational instabilities 
p9| and references therein). At low metallicity at the epoch of violent star formation, however. 



see 



they could have formed. The possibility of energetic GRBs from collapses of such massive stars was 
studied in [4^ with negative conclusion about their ability to produce an energetic GRB. But we 
note here that physical processes in such stars are still far from full understanding and potentially 
such stars could be GRB progenitors. The weakness of GRB980425 in a nearby galaxy can be a nat- 
ural consequence of smaller upper masses of stars in regions of violent star formation at the present 
epoch. 
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